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O^. Abstract 



We theoretically investigate excitation properties in the pseudogap regime of a trapped Fermi 
gas. Using a combined T-matrix theory with the local density approximation, we calculate strong- 
coupling corrections to single-particle local density of states (LDOS), as well as the single-particle 
local spectral weight (LSW). Starting from the superfluid phase transition temperature Tc, we clar- 
ify how the pseudogap structures in these quantities disappear with increasing the temperature. 
As in the case of a uniform Fermi gas, LDOS and LSW give different pseudogap temperatures T* 
and T** at which the pseudogap structures in these quantities completely disappear. Determin- 
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^^ \ ing T* and T** over the entire BCS (Bardeen-Cooper-Schrieffer)-BEC (Bose-Einstein condensate) 

crossover region, we identify the pseudogap regime in the phase diagram with respect to the tem- 
perature and the interaction strength. We also show that the so-called back-bending peak recently 
observed in the photoemission spectra by JILA group may be explained as an effect of pseudogap 
phenomenon in the trap center. Since strong pairing fluctuations, spatial inhomogeneity, and finite 
temperatures, are important keys in considering real cold Fermi gases, our results would be useful 
for clarifying normal state properties of this strongly interacting Fermi system. 

PACS numbers: 03.75.Ss,05.30.Fk,67.85.-d 



I. INTRODUCTION 



The pseudogap temperature is a fundamental quantity in strongly interacting Fermi su- 
perfluids, such as high-T^ cuprates[l|-l3| and superfluid ^°K and ^Li Fermi gases j#^17|. The 
normal state region below this characteristic temperature is referred to as the pseudogap 
regime, where anomalies are seen in various physical quantities, such as a gap-like struc- 
ture in single-particle density of states. As the origin of the pseudogap regime, the so-called 
preformed-pair scenario has been extensively discussed in high-Tc cuprates 



231. However, 



because of the complexity of this system, other scenarios have been also proposed, such as 
antiferromagnetic spin fluctuations 



24 



25|, and a hidden order J26|. In contrast, in cold 



Fermi gases, because of the 



simp 



icity of the system, the origin can be uniquely identified as 
strong pairing fluctuations |27N37l|. Thus, this system would be very useful for the study of 
the preformed-pair scenario. 

In considering the pseudogap phenomenon, we note that this is a normal state phe- 
nomenon, free from any phase transition at the pseudogap temperature. Because of this, 



the pseudogap temperature may depend on what we measure. In a previous paper 29|] for a 
uniform Fermi gas, we showed that the pseudogap temperature T* which is defined as the 
temperature at which the pseudogap structure disappears in the single-particle density of 
states is different from the pseudogap temperature T** determined from the single-particle 
spectral weight. While one finds T* > T** in the weak-coupling BCS regime, T** becomes 
higher than T*, as one passes through the BCS-BEC crossover region. 

Since a cold Fermi gas is always trapped in a harmonic potential, spatial inhomogeneity is 
also a key in considering the pseudogap problem of a real Fermi gas. Indeed, in the crossover 
region at Tc, it has been shown [30| that, while a clear pseudogap structure can be seen in 
the single-particle excitation spectrum in the trap center, one only sees a free-particle-like 
dispersion around the edge of the gas. Since the recent photoemission-type experiment 
developed by JILA group does not have spatial resolution, the observed spectra correspond 
to the sum of spatially inhomogeneous excitation spectra. Including this, we showed |30t] 
that the anomalous photoemission spectra observed at T^ can be theoretically reproduced. 
We briefly note that the importance of such inhomogeneous pseudogap phenomena in cold 
Fermi gases has been also pointed out in Ref . [34 1 . 

In this paper, we theoretically identify the pseudogap regime of a trapped Fermi gas 



in the BCS-BEC crossover region. Extending our previous work at Tc[30] to the region 
above Tc, we calculate the single-particle local density of states (LDOS), as well as the 
single-particle local spectral weight (LSW), within the framework of a combined strong- 
coupling T-matrix theory with the local density approximation (LDA). While the LDOS 
is suitable for the understanding of the meaning of 'pseudo'-gap, the LSW is closer to 
the photoemission spectrum. Starting from Tc, we show how the pseudogap structures 
in these quantities gradually disappear, as one increases the temperature. We define the 
pseudogap temperature T* as the temperature at which the pseudogap structure completely 
disappears in LDOS. In the same manner, we also define another pseudogap temperature 
T** for LSW. We determine both T* and T** in the entire BCS-BEC crossover region. Using 
these two pseudogap temperatures, we identify the pseudogap region in the phase diagram 
of cold Fermi gases. We also examine the photoemission spectrum, and explain the so-called 
back-bending behavior observed in ^"^K Fermi gases from the viewpoint of inhomogeneous 
pseudogap effect. 

This paper is organized as follow. In Sec. HIl we explain our formulation. Using the 
combined T-matrix theory with LDA, we first determine the chemical potential above Tc 
in the BCS-BEC crossover. We then calculate LDOS, as well as LSW, above Tc. Here, 
we also calculate the photoemission spectrum for a trapped Fermi gas within LDA. In 
Sec. Illlt we examine pseudogap effects on the LDOS and LSW, focusing on their temperature 
dependences. After determining the pseudogap temperatures T* and T**, we identify the 
pseudogap region in the phase diagram of trapped Fermi gases above Tc. In Sec. IV, we 
discuss the photoemission spectrum. We clarify how the pseudogap affects this quantity. 
Throughout this paper, we take h = ks = I- 

II. FORMULATION 

We consider a two-component Fermi gas, described by the BCS Hamiltonian, 

p,a P,p',q 

where c^^ is the creation operator of a Fermi atom with pseudospin a =t, |, which de- 
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scribe two atomic hyperfine states, ^p = £p — fJ^ = ^ ~ /^ is the kinetic energy, mea- 
sured from the chemical potential /x (where m is an atomic mass). Here, A^ is the to- 
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FIG. 1: (a) Self-energy $]p(ict;„), and (b) particle-particle scattering matrix rq(iz/„), in the T-matrix 
approximation. The solid and wavy lines represent the noninteracting Green's function Gp(ia;„) 
and the pairing interaction —U, respectively. 

tal number operator of Fermi atoms. —U{< 0) is an attractive interaction which can be 
tuned by a Feshbach resonance 38| . U is related to the s-wave scattering length a^ as 
4:TTas/m = —U/[l — U ^p'' l/(2ep)] (where Uc is a high-energy cutoff) [39]. As usual, we con- 
veniently measure the interaction strength in terms of the inverse scattering length (fc^a^)"^, 
where kp is the Fermi momentum. We will include effects of a trap later. 



We treat pairing fluctuations within the standard T-matrix theory |20l. |22| . |23| . |29|-|31|. 
Despite its simplicity, previous studies have shown that this strong-coupling theory gives 
quantitatively reliable results for pseudogap phenomena[30|. Within this framework, the 
single-particle thermal Green's function is given by 



Gp{iujn) 



Glitu^ 



(2) 



-Sp(za;„)' 

where Un is the fermion Matsubara frequency. Gp^iUn) = l/[i(jJn — ^p] is the single-particle 
Green's function for a free fermion. The self-energy part Tip{iuJn) involves effects of pairing 
fluctuations within the T-matrix approximation, which is diagrammatically given in Fig. [H 
Summing up these diagrams, we obtain 



T.p{iuJn) = T^ rq(iz/„)G'J_p(iz/„ - iun), (3) 

where z/„ is the boson Matsubara frequency. The particle-particle scattering matrix TgliUn) 



is given by 
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(4) 



The pair-correlation function 

P,^n 

1 - f{ip+q/2) - f{i~p+q/2) 



-lUJr, 



E 



ip-^q/2 + i-p+q/2 — il^n 



(5) 



describes fluctuations in the Cooper channel, where f{^) is the Fermi distribution function. 
We now include effects of a harmonic trap V{r) = mu^j.r'^/2 within LDA (where cutr 
is a trap frequency). This extension is simply achieved by replacing the Fermi chemical 
potential fi by the LDA expression fi{r) = yU — V{r) 40| in Eqs. (I2])-(|5]). In this paper, we 
explicitly write the variable r for LDA quantities. For example, the LDA single-particle 
Green's function is written as 

1 



Gp{iujn, r) 



i^n-^p{r) -Sp(zw„,r)' 



(6) 



where ^p{r) = Ep — fi{r). 

To examine the pseudogap phenomenon in a trapped Fermi gas, we consider the single- 
particle local density of states (LDOS) p{uj,r), as well as the single-particle local spectral 
weight (LSW) A{p,u,r). Their LDA expressions are given by 

p{u, r) = y^ \m[Gp{iUn -^ u + i6, r)], (7) 
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A{jp,u, r) = lm.[Gp{iuJn -^ u + i5,r)]. 



TT 



(8) 



The former can be written as p{u,r) = ^^pAlpjUjv), so that LSW may be viewed as the 
momentum resolved LDOS. 



In 



group 



:his paper, we also consider the photoemission-type experiment developed by JILA 



27 



28| . When the ^-spin state is coupled with another hyperfine state |3) (7^ | t)) I i)) 



by radio-frequency (rf) pulse [27 



the photoemission spectrum is given by 



p, Q) = '^l drA{p, Ur) - ^, r)f{Ur) - fi). 



(9) 



(We summarize the derivation of Eq. (Q in the Appendix.) Here, V = AirRp/S (where 
Rp = A/2/i/(r7ici;^j.) is the Thomas- Fermi radius [40(]), and tp is a couphng between the ^-spin 
state and |3). Since the current experiment has no spatial resolution, we have taken the 
spatial average in Eq. ([9]). Apart from this spatial integration, the photoemission spectrum 
is closely related to the spectral weight A{p,u,r). 



To directly compare our results with the experimental data 27, 
Eq. Q as 



we slightly modify 



A{p, Uj)f\uj) = /avc(P, fi ^ ^p - W) 






drA{p, u - fi{r),r)f{uj - n{r)). 



For a free Fermi gas, Eq. flTUl) is evaluated to give 

S7r3/2f2 

A(p,a;)/(u;) = -:;?UfeLi 



e-nSico-U 



(10) 



(11) 



2(/3/i)3/2^^3/2l 

where ljis{z) is the polylogarithm of z in the order s. The peak energy of Eq. (Ilip gives the 



single-particle dispersion u = C,p. At T = 0, Eq. f lTU]) reduces to 

3/2 

S{cu-Q9{-co) 



A{p,u)fiuj) = 2n4 



JILA's experiments 27|, 



(12) 



have also examined the occupied density of states, defined by 

p{uj)f{uj) = — ^ Y^ /ave(P, fi ^ ^p 



27rt^ 



-u 



(13) 



For a free Fermi gas, Eq. ( fT3l) gives 

3m3/2 
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4v^7r3/2 (/3/x)3/2 



Li3/2(-e"^"')v/wT7^0(w + /i). 



At T = 0, Eq. (dH) reduces to 



p(w)/(w) 
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Vw + /i^(u; + /i)^(-u;). 



(14) 



(15) 



Equation ( [T5|) is just the ordinary density of states oc ^/u^+~j2 multiplied by [wp/^ when 
a; < and u + ft >0. 

In order to calculate Eqs. ([7]), ([8]), (TTOj) . and (TT3|) . we need to determine the Fermi 
chemical potential /i from the equation for the total number N of Fermi atoms. The LDA 
number equation is given by 



A^ = 2T f drJ^Gpi 



(16) 
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FIG. 2: Calculated Fermi chemical potential /x above Tc. ep = (3A^)^'^ti;tr is the Fermi energy 
of a trapped Fermi gas, and the interaction is measured in terms of the inverse scattering length 
{kFas)~^, where kp is the Fermi momentum. We will use this result in calculating Eqs.([7])-([9]) in 
Sees. Ill and IV. The upper left and right panels show Tc and /i(T'c), respectively. 



The calculated chemical potential fi{T > T^) is shown in Fig. [2J We briefly note that 
the LDA superfluid phase transition temperature T^ is given as the temperature at which 
the Thouless criterion 






4l| is satisfied in the trap center (r = 0)[l6|. The resulting LDA Tc 



equation is given by 



r^=o(^i^n = o,r = o)-i = o. 



(17) 



III. PSEUDOGAP TEMPERATURES DETERMINED FROM LDOS AND LSW 



Figure [3] shows the LDOS p{uj,r) in the unitarity limit {{kpas)^^ = 0). At T^, panel (a) 
shows that LDOS in the trap center (r = 0) has a large dip structure around u = 0. Since 
the superfluid order parameter vanishes at Tc, this pseudogap structure purely arises from 



strong pairing fluctuations in the unitarity regime |29|-|31|. The pseudogap (dip) structure 



gradually disappears, as one goes away from the trap center, because pairing fluctuations 





FIG. 3: (Color online) Calculated local density of states (LDOS) p{uj,r) above Tc in the unitarity 
limit [{kpas)^^ = 0). The fine peaks seen in LDOS at large r are a computational problem in the 
momentum integration of A{p,uj,r), involving (5-functional peaks as A{p,uj,r) ~ 6{oj — ^p{r)). 

become weak around the edge of the gas (r ~ Rp) due to the low particle density. 

With regard to inhomogeneous pairing fluctuations, we note that pairing fluctuations are 
described by the analytic continued particle-particle scattering matrix, 

Ana^ 1 



^q{i^n — )■ w + i5, r) 



m 



1 + 



A-Kaa 



Ilq{iVn -^^ + i5,r)-Y.p^ 



(18) 



where 



^q{i^n — )■ W + «5, r) = 2_\ 



l-/(ep+./2(r))-/(^p+,/2(r)) 



ip+q/2{r) +i-p+q/2{ 



{u + i5) ' 



(19) 



Because of the Thouless criterion in Eq. 01 7p . in the trap center (r = 0), Eq. (1181) with 
w = qr = diverges at T^. On the other hand, because the LDA chemical potential ^{r) = fx— 
V{r) = fi — 'muj'^r'^/2 decreases as one goes away from the trap center (r > 0), 11^=0(0, r ^ 0) 
becomes small. Thus, low energy and low momentum pairing fluctuations described by 
r<j~o(i^ ~ 0,r > 0) also become weaker for larger r. In particular, around the edge of the 
gas cloud, panel (a) shows that the LDOS becomes close to that of a noninteracting Fermi 
gas. 



p{u},r) 



3 

V2n^ 
8 



^UJ + fi{r). 



(20) 



(k^a^'-l 



fi 2 



Q^^y^ 



<k,a;'=l 



25 




0.25 



FIG. 4: (Color online) Calculated intensity of local spectral weight (LSW) A(p, u) in the trap 
center r = 0. The intensity is normalized by the inverse Fermi energy ej, . 

As one increases the temperature, the pseudogap structure in LDOS gradually disappears 
from the outer region of the gas cloud, as shown in Figs. |3](b)-(d). Defining the pseudogap 
temperature T* as the temperature at which the dip structure in the trap center disappears, 
one obtains T* = l.lTc in the case of Fig. [3l Above the pseudogap temperature T*, Figj3l^d) 
shows that the overall structure of LDOS becomes close to that for a free Fermi gas in Eq. 
( 120|) . although, in the trap center, LDOS around the threshold energy (w = —fi{r)) is found 
to be still affected by pairing fluctuations to deviate from p{u, r) ~ a/w + /^(r). 

Figure m shows the LSW A{p,u,r) at r = 0. In the BCS side (panels (a)-(c)), although 
a pseudogap structure can be slightly seen at p ~ /cf in panel (a), it soon disappears at 
higher temperatures. (See panels (b) and (c).) As discussed in Ref. J30l], the pseudogap 
effects on LSW soon disappears, as one moves away from the trap center (although we do 
not explicitly show the spatial dependence of LSW in the BCS side in this paper). Thus, 
Fig. Ill^a)-(c) indicate that LSW is not useful for the observation of pseudogap phenomenon 
in the BCS side. 

In the unitarity limit, we see a clear pseudogap structure in LSW at r = 0, as shown in 



Figs. |l]^d)-(f). In particular, at Tc, although the superfluid order parameter vanishes, the 
overall pseudogap structure is similar to the superfluid gap structure in the BCS spectral 
weight, given by, 

A{p, u,r = 0)= ul6{u - Ep) + vl6{uj + Ep), (21) 

where Ep = ^^^ + A^ Up = ^/{l + ^p/Ep)/2, and Vp = ^{l - ip/ Ep)/2. Here, A is the 
BCS superfluid order parameter. In the BCS case, the double peak structure at a; = ±£'p 
may be regarded as a result of a coupling between particle and hole excitations by the 
superfluid order parameter A. Indeed, the diagonal component of the BCS Green's function 
can be written in the form 



Gu{p,iUJn) 



c^^ + e^ + A^ 

^ - (22) 



[zUn - ^p 



i^n + ip 

Noting that l/(^c^;„ — ^p) and l/{iuJn +^p) represent the particle and hole Green's functions, 
respectively, A^ in Eq. (1221) is found to work as a coupling between the particle branch and 
the hole branch. In the pseudogap case at Tc, since the particle-particle scattering matrix 
^q{iT-'n) diverges at z/^ = q = (Thouless criterion), one may approximate the self-energy in 
Eq. ([3D to 

T.p{lUJn) ^TJ2 rq(^'^n) X G'Lpi-lUJn). (23) 

When one substitutes Eq. (123!) into Eq. ([2]), the resulting expression has the same form as 



Eq. fl22|) where A^ is replaced by the pseudogap parameter A'^^ = ^TJ2qu„ r^(zz/„)|23l. |29|. 
That is, pairing fluctuations described by Tg^iUn) induces a particle-hole coupling, leading to 
the pseudogap structure seen in Fig. Hl^d). Since pairing fluctuations become weak at higher 
temperatures, the double peak structure in the spectrum becomes obscure to eventually 
disappear, as shown in Figs. |l](e) and (f). 

Besides the particle-hole coupling, pairing fluctuations also lead to flnite lifetime of quasi- 
particle excitations. Because of this, the LSW in the unitarity limit exhibits broader spectra 
compared with that in the BCS regime, as shown in Figs. H] (a)-(f). 

In the strong-coupling BEG limit, the system reduces to a gas of two-body bound 
molecules, so that single-particle excitations are simply described by dissociations of these 
two-body bound states. Then, noting that the single-particle spectral weight in the negative 
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FIG. 5: (Color online) Calculated intensity of LSW A{p,uj,r) in the unitarity limit {{kpas) ^ = 0). 
The intensity is normalized by the inverse Fermi energy ej, . 



energy region {u < 0) physically describes hole excitations, we expect that the intensity of 
the lower spectral branch becomes weak in this two-body regime. Indeed, the lower peak is 
found to be very weak in Fig. |4l^g)-(i). We briefly note that the peak width of the upper 
branch in the BEC regime is sharper than the case of unitarity limit shown in Figs. Hl^d)- 
(f), which is due to the fact that this branch simply describes the dissociation energy of a 
two-body bound state in the BEC limit. 

Figures E] and El respectively, show the LSW in the unitarity limit {{kpas)^^ = 0) and 
in the BEC regime ((fci^a^)"^ = 1). We find that the pseudogap effect is not so remarkable 
around the edge of the gas cloud. We also find that the r-dependence of LSW in Fig. [5] 
is similar to the temperature dependence of this quantity shown in Figs. H] (d)-(f). To 
understand this similarity, we recall that the spatial dependence of LSW is dominated by 
the LDA chemical potential fi{r) = /i — V{r), which decreases as r increases. The chemical 
potential /i also becomes small with increasing the temperature, as seen in Fig. O As a 
result, the increase of r and the increase of the temperature lead to the similar effect on 
LSW seen in Figs. Hd)-(f) and Fig. 15142 1. 
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FIG. 6: (Color online) Calculated intensity of LSW A{p,uj,r) in the BEC side {{kpasy^ = 1). 
The intensity is normalized by the inverse Fermi energy ej, . 

The lower peak in the BCS spectral weight in Eq. (^T^ at a; = —Ep is a downward curve in 
the BEC regime where // < 0. Such a behavior is also seen in the lower peak line, as shown in 
Fig. m^g). However, Figs. 111(h) and (i) show that this downward dispersion gradually changes 
into an upward one with increasing the temperature. The resulting upward dispersion of 
the lower branch is similar to the hole branch in the presence of a Fermi surface {^ > 0), so 
that this phenomenon would have to do with the hole-type character of the lower branch. 
At high temperatures, since unpaired fermions are thermally excited to occupy the upper 
branch of the excitation spectrum, these free-particle-like atoms are expected to give the 
upward spectrum in LSW. 

As in the unitarity limit, the r-dependence of LSW in Fig. [6]is similar to the temperature 
dependence of this quantity in Figs. H] (g)-(i). However, in contrast to the behavior in the 
unitarity limit, the upper peak becomes slightly broad as increasing the temperature. This is 
clearly seen at the edge of the trap (r = Rp) in Figs, [nid)-(f). The upper peak becomes broad 
as the upward dispersion of the lower branch becomes remarkable, so that this broadening 
would be also caused by thermally excited fermions. Since the effective temperature increases 
as one moves away from the trap center, this effect is most remarkable at the edge of the 
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FIG. 7: (Color online) Two pseudogap temperatures T* and T** determined from LDOS and LSW, 
respectively. The solid line shows the superfluid transition temperature Tc, below which the system 
is in the superfluid phase (SF). In this figure, we also plot 2|^(Tc)| in the normal state of the BEC 
regime where //(Tc) < 0. Since 2\fi\ equals the dissociation energy of a two-body bound state in 
the BEC limit, T ~ 2|^(Tc)| physically gives the characteristic temperature where two-body bound 
molecules are formed. Thus, the right side of 2|^(Tc)| may be viewed as a normal-state molecular 
Bose gas (NB), rather than a Fermi gas. The pseudogap regime (PG) is the region surrounded 
by Tc, 2|/i|, and T* or T**. The left side of the pseudogap regime is the normal Fermi gas regime 
(NF), where strong-coupling effects are not crucial. We note that T*, T**, and 2^(Tc), are all 
characteristic temperatures, free from any phase transition. 



trap. 

When we define another pseudogap temperature T** as the temperature at which the 
double peak structure in LSW completely disappears, it does not coincide with the pseudo- 
gap temperature T* determined from LDOS, as shown in Fig. [71 As in the homogeneous 
case|29|, while T* > T** in the BCS side, T** becomes higher than T* when (fepa,)"^ > -0.1. 
Since the pseudogap is a crossover phenomenon, the pseudogap temperature may depend 
on what we measure. 

In Fig. [TJ we also plot 2|/i(Tc)| in the BEC regime where /i(Tc) < 0. Since 2|/i| in the 
BEC limit equals the binding energy of a two-body bound state, this line physically gives the 
characteristic temperature where two-body bound molecules starts appearing, overwhelming 
thermal dissociation. Thus, the right side of T = 2\fi{Tc)\ in Fig. [7]may be regarded as a gas 
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of two-body bound molecules, rather than a Fermi atom gas. Including this, we identify the 
pseudogap regime as the region surrounded by Tc, 2|/i(Tc)|, and the pseudogap temperature 
T* or T** (which depends on which we measure, LDOS or LSW). 

Although the phase diagram in Fig. [7] is very similar to that for a homogeneous Fermi 
gas 29| , we note that the pseudogap regime in the trapped case is narrower than that in 
the homogeneous case, in the sense that the ratios T* /Tc and T**/Tc in the former case 
are smaller than those in the latter. This is because, while pairing fluctuations are strong 
everywhere near Tc in the uniform case, strong pairing fluctuations are restricted to the 
spatial region around the trap center in the trapped case. As a result, the pseudogap 
phenomenon is somehow weakened by the outer region of the gas cloud. 

We find in Fig. [7]that the pseudogap region is very narrow in the BCS side {{kpag)^^ < — 
0.1), indicating difficulty of observing the pseudogap in this weak-coupling regime. On the 
other hand, in the crossover region ((fcpas)"^ > —0.1), although the pseudogap temperature 
T* is still not so high, T** is found to give a large pseudogap regime. Since the photoemission 
spectrum /ave(9, ^) in Eq- & is close to the spectrum weight A{p, u, r), the large pseudogap 
region in the BEG side makes us expect that the photoemission-type experiment would be 
useful in observing the pseudogap in trapped Fermi gases. In Sec. IV, we will check this 
expectation by explicitly calculating the photoemission spectrum in the BEG side. 

To see inhomogeneous pairing fluctuations in a simple manner, it is convenient to 
divide the atomic density profile n{r) into the sum of the free fermion part n-pir) = 
2r^ ^^ Gp{iUn,r)e'^'^"^ and the fiuctuation correction 2nB{r), where 

nB(r) = tJ2 [Gp{^cOn, r) - ^^(za;,, r)] e'^-' . (24) 

In the strong-coupling BEG regime where /i ^ —ep, n-sir) reduces to the number of tightly 



bound molecules ll| . As shown in Fig. [HI the fiuctuation contribution n-Q{r) always takes 
a maximum value in the trap center, as expected. While the density profile n{r) is always 
dominated by np{r) above Tc in the BGS regime (panel (a)), n-sir) in the BGS-BEG crossover 
region is remarkably enhanced around the trap center near Tc, as shown in panels (b)-(d). 

When (/cpfls)"^ = 1 shown in panel (d), the density profile is dominated by nB{r). In 
addition, we also see a cusp in nB(r) around r = at Tc. Since this cusp structure is 
characteristic of the LDA density profile of a Bose gas at Tc[40], we may regard the superfiuid 
phase transition in this regime as a BEG of molecular bosons described by nB(r). 
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FIG. 8: (Color online) Density profile of free fermion component nF(r) and fluctuation correction 
nB(r) in the BCS-BEC crossover. 

IV. PHOTOEMISSION SPECTRUM AND PSEUDOGAP EFFECTS 

Figure |9] shows the photoeinission spectrum above Tc. (We also show in Fig. [10] the same 
plots as in Figs. |9t^a)-(c) in the logarithmic scale, in order to clearly show the peak structure.) 
Although we do not explicitly show results in the weak-coupling BCS side ((/cfOs)"^ < 0), 
as expected from the phase diagram in Fig. [TJ the photoemission spectrum in this regime 
is essentially the same as that for a free Fermi gas. That is, one only sees a sharp peak line 
along the free particle dispersion uj + fi = p^ /{2m). On the other hand, in the crossover 
and BEG regime shown in Fig. |9l the photoemission spectrum at Tc has an upper sharp 
peak at w + /i ~ p"^ /(2m) (upper dashed line) and a lower broad peak (lower dashed line). 
As discussed in Ref. |30|, this double peak structure originates from the inhomogeneous 
pseudogap phenomenon seen in LSW. That is, the upper sharp peak dominantly arises from 
the LSW around the edge of the gas cloud where pairing fluctuations are weak, so that the 
upper peak position is close to the free particle dispersion. On the other hand, the lower 
broad peak directly reflects the lower branch in the pseudogapped LSW in the trap center 
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FIG. 9: (Color online) Calculated intensity of photoemission spectrum p'^A{p,Lo)f{Lo) above Tc. 
The intensity is normalized by 27rtp/(2?7i). The spectrum has a sharp peak along the free particle 
dispersion oj = S,p (upper dashed line). In addition, it also has a lower peak line, which is shown 
as the lower dashed line. Since the intensity of the lower peak is much weaker than the upper 
one, in some panels, one cannot see the lower one in the intensity plot. Solid circles represent the 
experimental data observed in Ref. [27|. 

where pairing fluctuations become strong near T^.. 

As discussed in Sec. Ill, the appearance of the lower branch in LSW is a direct conse- 
quence of the particle-hole coupling induced by pairing fluctuations. Because of the similarity 
between this coupling effect and that induced by the superfluid order parameter below T^ 
the lower peak line seen in Fig. [9](a) exhibits a back-bending behavior, being similar to 
the BCS hole excitation spectrum — ^y^^y^ (although the momentum at which the back- 



bending occurs is different from p = y/2mj2 < kp in the present case). In the BCS state, the 
back-bending behavior of the hole excitations spectrum — ^J^p + A^ vanishes in the strong- 
coupling BEC regime where /i < 0. This tendency can be also seen in the photoemission 
spectrum at Tc, as shown in Figs. M^a), (d), and (g). In Figs. |9t^a) and (g), we find that the 
lower peak lines evaluated in the unitarity limit and BEC regime agree well with the recent 
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FIG. 10: (Color online) Same plots as in Figs. [9)^a)-(c), in the logarithmic scale, in Ref. 

experiments on "^^K Fermi gases[27|. This indicates that the observed back-bending behavior 
may be understood as a signature of the pseudogap effect. Furthermore, these agreements 
demonstrate the vahdity of our strong-couphng theory for the pseudogap physics in cold 
Fermi gases. 

Since the pseudogap temperature T** is close to Tc in the unitarity limit (See Fig. [71), the 
lower spectral peak line in the low momentum region p/kp < 1.5 soon vanishes with increas- 
ing the temperature, as shown in Figs. [9]^b) and (c). However, even at T = 1.98Tc, panel 
(c) shows that the lower peak line still remains in the high momentum region p/kp > 1.5, 
although the peak is very broad and the intensity of the peak is weak, as shown in Figs. 
[TOT b) and (c). Since this temperature is much higher than T** = 1.26Tc, this lower peak line 
remaining at high momenta would be nothing to do with the pseudogap effect. Instead, as 
pointed out in Ref. 43l-l45| . this lower spectral peak is considered to arise from the universal 
behavior of a Fermi gas with a contact interaction. In this regard, we recall that, below 
T**, the back-bending behavior of the lower peak is due to the pseudogap effect, originating 
from the double-peak structure of the LSW in the trap center. On the other hand, such a 
double-peak structure of the LSW is absent above T**, so that the back-bending behavior in 
the high temperature regime is purely due to the contact interaction. That is, the origin of 
the back-bending behavior of the lower peak line changes around the pseudogap temperature 

When (/cpfls)^^ = 0.4 shown in Figs. [9]^d)-(f), the lower peak line with the back-bending 
behavior can be clearly seen even at = 1.5Tc, because of the wide pseudogap regime in this 
case. (See Fig. [71) As one increases the temperature, the bending position moves to higher 
momenta, which is consistent with the recent experiment above T^. 28j. 
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FIG. 11: (Color online) Calculated photoeniission spectrum p'^(A(p,u})f(uj)) affected by finite 



271 . (a) {kpas 



0. (b) 



resolution, given by Eq. (j26p . We take T = Tc, and E^- = 0.2Ep 

(/cfOs)"^ = 0.4. (c) {kpUs)"^ = 1. Panel (d) shows the effect of finite resolution on the occupied 

density of state, given by {p{uj) f {lj)) = ^ (A(p, u;)/(a;))/(27rtp). In this figure, experimental 



data 



27t | are shown as solid circles. In panel (d), we have offset the curves by 0.1. 



We briefly note that, although the region at (fepOs)"^ = 1 may be regarded as a molecular 
Bose gas in our phase diagram in Fig. [TJ we still see the lower branch in Figs. [9](g)-(i). This 
indicates that many-body effects still contribute to the pair-formation to some extent even 
when {kpas)~^ = 1. As discussed in Sec. Ill, the upward curve of the lower peak line in 
the LSW reflects the hole-like character of this branch. This behavior is also seen in Figs. 
|9](g)-(i), where the lower downward peak line around /c ~ gradually becomes an upward 
one with increasing the temperature. 



Since the observed spectra are always affected by a finite energy resolution [27|, |28|, it is 
interesting to see how our results are modified by this smearing effect. To briefly examine 
this, we consider 



{A{p,u)f{u)) 



'27: E, 



dz A{p,z)f{z)e 



-{z-uj)yi2E?) 



(25) 
(26) 
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FIG. 12: (Color online) Calculated occupied density 



of states p{uj)f{uj) above Tc. In panels (a) 



271]. The magnitudes of the experimental data 



and (c), solid circles are experimental data in Ref. 

are normalized so that their peak heights can coincide with the peak heights of our results at Tc. 

We have offset the curves by 0.1 in panels (a) and (b), and by 0.05 in panel (c). 

where E^ describes a finite energy resolution. We then find in Fig. [11] that, although the 
resulting spectral peaks are broadened by the finite resolution, the overall spectral structure 
remain unchanged. Thus, the finite energy resolution {Er = 0.2£:f) is not a serious problem 
for the observation of the pseudogap phenomenon in the BCS-BEC regime of a cold Fermi 

gas. 



Figure [T2] shows the occupied density of states p{uj)f{uj) in Eq. (1T3|) . In the unitarity limit 
(panel (a)), a broad peak only appears above Tc, so that we cannot see a signature of the 
pseudogap. However, in the BEC side, panel (b) clearly exhibits a double peak structure 
above Tc, reflecting the pseudogap in single-particle excitations. This structure becomes 
more remarkable in panel (c). 

We compare the calculated occupied density of states with the experiment results on a 
^°K Fermi gasJ27|. In the unitarity limit, the single peak structure in the unitarity limit well 
agrees with the observed occupied density of states, as shown in Fig. [TW a). In the BEC 
regime (panel (c)), although the relative peak height between the upper and lower peak in 
our result is opposite to the experimental result, their peak positions are in good agreement 
with the experiment. These agreements still hold, even when one takes into account the 
experimental finite energy resolution, as shown in Fig. ITTT d). Thus, the present T-matrix 
theory is found to well describe strong-coupling effects on the occupied density of states. 
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V. SUMMARY 

To summarize, we have discussed the pseudogap phenomenon and effects of a harmonic 
trap in the BCS-BEC crossover regime of an ultracold Fermi gas. Extending our previous 
work at T^ to the region above T^, we have calculated single-particle local density of states 
(LDOS), as well as the single-particle local spectral weight (LSW), within the framework of 
the combined T-matrix theory with the local density approximation. We clarified how the 
pseudogap structures in these quantities gradually disappear with increasing the tempera- 
ture. From these, we introduced two pseudogap temperatures T* (for the density of states) 
and T** (for the spectral weight). 

As in the case of a uniform Fermi gas, T* does not coincide with T**. While T* > T** 
is obtained in the weak-coupling BCS side, one finds T* < T** in the strong-coupling EEC 
side. This means that the pseudogap temperature depends on what we measure. 

We have also examined the photoemission spectrum in the pseudogap regime. At present, 
the photoemission-type experiment has no spatial resolution, so that the observed photoe- 
mission spectra are spatially averaged ones. Although the spatial average smears inho- 
mogeneous pseudogap effects to some extent, we showed that the double peak structure, 
which is characteristic of the pseudogap effect, still survives even when such a smearing 
effect is taken into account. The photoemission spectra, as well as the occupied density of 
states, calculated in the crossover and BEC regime agree well with the recent experiments 
on *^°K Fermi gases, so that the combined T-matrix theory with LDA used in this paper is 
found to be a powerful theory to study the pseudogap physics in cold Fermi gases. Since 
the pseudogap temperature, which we have determined in this paper, is a key quantity in 
the pseudogap physics, the experimental determination of this characteristic temperature 
would be useful for the further understanding of strong-pairing fluctuations existing in the 
BCS-BEC crossover regime of a cold Fermi gases. 
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Appendix A: Derivation of Eq. ([9]) 

In this appendix, we explain the outhne of the derivation of Eq. (|9]) . In the photoemission- 
type experiment [27i], atoms in one of the two hyperfine states (= | t)) are transferred to 
an unoccupied hyperfine state |3) (7^ | t)i I i)) by rf-pulse. Although this experimental 
procedure is essentially the same as the rf-tunneling current spectroscopy, one can safely 
ignore the so-called final state interaction in the recent photoemission-type experiment on 
^°K 27|. Noting this, we consider the model Hamiltonian H = H + H^ + Ht, where H is 
given in Eq. ([1]), and 

i^3-/i3A^3 = 5^(£p + t^3-Ai3)&J,&p (Al) 

p 

describes the final state |3). Here, bp is an annihilation operator of a Fermi atom in |3), and 

U3 is the energy difference between | f) and |3). fi^ and N^^ are the chemical potential and 
the total number operator of atoms in the final stat e |3 ), respectively. The transition from 



I t) fo |3) is described by the tunneling Hamiltonian [l6l. l46l. |47| 

HT = tFY^ (e'*""*4+.,Cfct + H.c.) , (A2) 

fc 

where tp is a transfer matrix element between | f) and |3). Ql and Ul represent the 

momentum and energy of the rf-pulse, respectively. As usual, we first consider a uniform 

Fermi gas, and then include effects of a trap by replacing fi^ with the LDA expression 

fJ'sir) = /is - V{r). 

The photoemission spectrum is conveniently described by the rf-tunneling current I from 

the initial state | t) to the final state |3), given by 

I = {N^)=t{[H + H3 + HT,N,]) = {J), 

J = _,t^ ^ [e— -*6i^^^c,t - H.c] . (A3) 

k 

Here, N^ is the total number operator of Fermi atoms in the final state |3). Within the 



linear response theory 48|] in terms of the tunneling Hamiltonian Ht in Eq. flA2p . Eq. flA3P 

reduces to 

I = -t f dt' {[J{t),HT{t')])e''', (A4) 

J —00 
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where Hrit) = e^(^+^3)i^^e-^(^+^3)*, and J{t) = e^(^+^3)ije-*(^+^3)t. 

Equation flA4p can be conveniently evaluated from the corresponding thermal Green's 
function by analytic continuation. That is, when one introduces the correlation function in 
the Matsubara formalism, 

A(il^n) = ^ X^ Gk^{iUn)Gk+qr^,3{i(^n + il^n) , (A5) 

k,LU„ 

(where Gp^^^iUn) = '^/[ioon — {Sp — ^3)] is the single-particle Green's function for the final 
state 1 3)), Eq. flA4p is given by 

I = -2 lm[A{iiyn ^ ^ + H - Us + id)]. (A6) 

Here, ^2 = wl — ws is the rf-detuning. Carrying out the summation with respect to the 
Matsubara frequencies in Eq. (lASp . one finds 

i{p,n) = 2ntlA{p,^p - n)f{^p - n). (at) 

In obtaining Eq. (]A7p . we have assumed that (1) the photon momentum is negligibly small 
{ql = 0), and (2) the final state |3) is initially empty {f{sp — /is) = 0). We find from Eq. 
(]A7p that the photoemission spectrum is related to the occupied spectral weight as 

I{p,Q^^p-u) = 27rtlA{p,uj)f{uj). (A8) 



At T = 0, since the Fermi distribution function f{uj) reduces to the step function, Eq. (lASp 
gives the ordinary spectral weight A{p, u) below 00 = 0. At finite temperatures, on the other 
hand, thermally excited quasiparticles contribute to the spectrum, so that the photoemission 
spectrum J(p, fi — )■ ^p — w) in Eq. (jASP has finite intensity even for cj > 0. 

The extension of Eq. flA7p to a trapped gas is achieved by replacing /i with the LDA 
expression ^{r) = /i — V{r), leading to the local photoemission spectrum, given by 

lip, n, r) = 27itlA{p, ep(r) - fi, r)f{^pir) - fi), (A9) 



where Cp{r) = Sp — fi{r). In the current experiment [271] , since the rf-pulse is always applied 
to the whole gas cloud, the observed spectrum involves contributions from all the spatial 
regions of the gas cloud. To include this, we take the spatial average of Eq. flA9p as 

/ave(p,fi) = ^ J drAip,Ur)-^,r)f{Ur)-n). (AlO) 
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[muD 



Here, V = 47ri?|./3 is a characteristic volume of the gas cloud, where Rp = A/2/i/(r 
is the Thomas- Fermi radius J40|. We emphasize that Eq. flAlOp gives a proper definition of 
the spatially averaged photoemission spectrum, being comparable to the observed spectra 
by JILA group. 



[1] 
[2] 

[3] 
[4] 



[5] 
[6] 

[7] 

[8] 

[9 

[lo; 

[11 

[12 
[13 
[14 

[15 
[16 
[17 
[18 



A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 75, 473 (2003). 

0. Fischer, M. Kugler, I. Maggio-Aprile, C. Berhod, and C. Renner, Rev. Mod. Phys. 79, 353 
(2007). 

P. Lee, N. Nagaosa, and X. Wen, Rev. Mod. Phys. 78, 17 (2006). 

For review, see Q. J. Chen, J. Stajic, S. N. Tan, and K. Levin, Phys. Rep. 412, 1 (2005); 

1. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008); S. Giorgini, L. 
Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008). 

C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 040403 (2004). 

M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Raupach, A. J. Kerman, and W. 

Ketterle, Phys. Rev. Lett 92, 120403 (2004). 

J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and J. E. Thomas, Phys. Rev. Lett. 92, 

150402 (2004). 

M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochini, C. Chin, J. Hecker Denschlag, and R. 

Grimm, Phys. Rev. Lett 92, 203201 (2004). 

D. M. Eagles, Phys. Rev. 186, 456 (1969). 

A. J. Leggett, Modern Trends in the Theory of Condensed Matter (Springer, Berlin, 1980). 

P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 195 (1985). 

C. A. R. Sa de Melo, M. Randeria, and R. Engelbrecht, Phys. Rev. Lett. 71, 3202 (1993). 

E. Timmermans, K. Furuya, P. W. Milonni, and A. K. Kerman, Phys. Lett. A 285, 228 (2001). 
M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo, and R. Walser, Phys. Rev. Lett 87, 
120406 (2001). 

Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89, 130402 (2002). 

Y. Ohashi and A. Griffin, Phys. Rev. A 67, 033603 (2003). 

Y. Ohashi and A. Griffin, Phys. Rev. A 72, 013601 (2005). 

M. Randeria, N. Trivedi, A. Moreo, and R. T. Scaletter, Phys. Rev. Lett. 69, 2001 (1992). 

23 



[19] J. M. Singer, M. H. Pedersen, T. Schneider, H. Beck, and H.-G. Matuttis, Phys. Rev. B 54, 

1286 (1996). 
[20] B. Janko, J. Maly, and K. Levin, Phys. Rev. B 56, R11407 (1997). 
[21] Y. Yanase and K. Yamada, J. Phys. Soc. Jpn. 70, 1659 (2001). 
[22] D. Rohe and W. Metzner, Phys. Rev. B 63, 224509 (2001). 

[23] A. Perah, P. Pieri, G. C. Strinati, and C. Castellani, Phys. Rev. B 66, 024510 (2002). 
[24] D. Pines, Z. Phys. B 103, 129 (1997) and references are therin. 
[25] A. Kampf and J. R. Schrieffer, Phys. Rev. B 41, 6399 (1990). 

[26] S. Chakravarty, R. B. Laughhn, D. K. Morr, and C. Nayak, Phys. Rev. B 63, 094503 (2001). 
[27] J. T. Stewart, J. P. Gaebler, and D. S. Jin, Nature 454, 744 (2008). 
[28] J. P. Gaebler, J. T. Stewart, T. E. Drake, D. S. Jin, A. Perah, P. Pieri, and G. C. Strinati, 

Nature Physics, 6, 569 (2010). 
[29] S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 80, 033613 (2009). 
[30] S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 82, 033629 (2010). 
[31] R. Watanabe, S. Tsuchiya, and Y. Ohashi, Phys. Rev. A 82, 043630 (2010). 
[32] Q. Chen and K. Levin, Phys. Rev. Lett. 102, 190402 (2009); C. Chien, H. Guo, Y. He, and 

K. Levin, Phys. Rev. A 81, 023622 (2010). 
[33] P. Magierski, G. Wlazlowski, A. Bulgac, and J. E. Drut, Phys. Rev. Lett. 103, 210403 (2009); 

P. Magierski, G. Wlazlowski, and A. Bukac. [a?Xiv:1103.438 2 (2011). 
[34] H. Hu, X.-J. Liu, P. D. Drummond, and H. Dong, Phys. Rev. Lett. 104, 240407 (2010). 
[35] A. Perali, F. Palestini, P. Pieri, G. C. Strinati, J. T. Stewart, J. P. Gaebler, T. E. Drake, and 

D. S. Jin, Phys. Rev. Lett. 106, 060402 (2011). 
[36] S. Su, D. E. Sheehy, J. Moreno, and M. Jarreh, Phys. Rev. A 81, 051604(R) (2010). 
[37] E. J. Mueller, Phys. Rev. A 83, 053623 (2011). 

[38] For review, see C. Chin, R. Grimm, P. Juhenne, E. Tiesinga, Rev. Mod. Phys. 82, 1225 (2010). 
[39] M. Randeria, in Bose-Einstein Condensation, edited by A. Griffin, D. W. Snoke, and S. 

Stringari (Cambridge University Press, New York, 1995), p. 355. 
[40] C. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Cases (Cambridge University 

Press, Cambridge, United Kingdom, 2002). 
[41] D. J. Thouless, Ann. Phys. (N.Y.) 10, 553 (1960). 
[42] The self-energy correction in Eq. ^ also has a temperature dependence which is independent 

24 



of the temperature dependence of /i. However, this temperature dependence is actually weak, 

so that this effect is not clearly seen in Figs. U] and [5l 
[43] S. Tan, Ann. Phys. (N.Y.) 323, 2971 (2008); S. Tan, ibid. 323, 2987 (2008); S. Tan, ibid. 323, 

2952 (2008). 
[44] W. Schneider, and M. Randeria, Phys. Rev. A (2010). 

[45] J. T. Stewart, J. P. Gaebler, T. E. Drake, and D. S. Jin, Phys. Rev. Lett. 104, 235301 (2010). 
[46] P. Torma and P. Zoller, Phys. Rev. Lett. 85, 487 (2000). 
[47] Y. He, Q. Chen, and K. Levin, Phys. Rev. A 72, 011602(R) (2005). 
[48] G. Mahan, Many Particle Physics (Plenum, New York, 1981). 



25 



